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Abstract 

Motivated by the recent works of us [1], we establish a correspondence be- 
tween /(i?)-gravity with the power-law entropy corrected holographic dark energy 
(PLECHDE) and new agegraphic dark energy (PLECNADE) models. We recon- 
struct corresponding /(i?)-gravities and obtain the equation of state parameters. 
We show that the selected /(i?)-gravity models can accommodate the transition 
from the quintessence state to the phantom regime as indicated by the recent ob- 
servations. 
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1 The theory of f(R) modified gravity 



The Friedmann equations in /(i?)-gravity are given by [2] 
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H 2 = p m + PR, 

(2H + 3H 2 ) = -(p m +p R ), 



where 
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with 



--f(R) + 3{H + H 2 )f'{R) - 18{4H 2 H + HH)f'(R) 
y {R) - { H + 3H 2 )f'(R)+ 

6(8H 2 H + 6HH + AH 2 + H)f"(R) + 3Q(HH + AH 2 H) 2 f"'(R) 



R = 6(H + 2H 2 ). 
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Here, p rn and p m are the energy density and pressure of the matter, respectively. Also p R 
and pr are due to the contribution of /(-R)-gravity. 

The equation of state (EoS) parameter of the curvature contribution is defined as [3] 



PR , . Pr + Pr 
u r = — = -H . 

PR PR 



(6) 



For a given a = a(t), by the help of Eqs. (3) and (4) one can reconstruct the /(i?)-gravity 
according to any DE model given by the EoS pn = Pr(pr) or p R = Pr(o) [4]. 
Here we consider the two classes of scale factors, the first one is given by [4, 5] 



a(t) = a (t s -ty h , t<t s , h>0. 
Using Eqs. (5) and (7) one can obtain 
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For the second class of scale factors defined as [4] 

a(t) = a t h , h>0, 

one can get 



, H = H 2 /h. 



H = 



h 



h 



-R 



1/2 



, H = -H 2 /h. 



(7) 
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_Q(2h - 1)' 

In sections 2 to 4 using the two classes of scale factors (7) and (9), we reconstruct different 
/(^-gravities according to the PLECHDE and PLECNADE models. 
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2 PLECH /(fl)-gravity model 

The holographic dark energy (HDE) density is given by [6] 



PA = 



3M 2 c 2 
L 2 : 



where c and M~ 2 = 8nG are the numerical constant and reduced Planck mass. Also 
L is the IR-cutoff of the universe. This model can be modified due to the power-law 
corrections to the entropy which appears in dealing with the entanglment of quantum 
fields in and out the horizon [7] . The power-law corrected entropy takes the form [7] 



b ~ AG 



1 - K n A 1 ~ a/2 



(12) 



where a is a dimensionless constant whose value is currently under debate and 

U \a/2-l 

K ° - jhhf-- (13) 

Also r c is the crossover scale. Using (12), the PLECHDE density was introduced as [8] 

-3c 2 f3 



(14) 



where in the absence of correction terms (a = (3 = 0) yields the well-known HDE density 
[7]- 

For a flat universe we have L = R h , where Rh is the event horizon defined as 



da 



(15) 



r°° at r c 
Rh = a — = a / -— . 

For the first class of scale factors (7) and using Eq. (8), the future event horizon Rh yields 

1 



Rh = a — 

Jt a 



*• dt t s -t 



h+1 h+l 



l 6h(2h+l) 
R 



(16) 



Here due to having R h > then —l<h. 
Inserting Eq. (16) into (14) one can get 



Pa = M 



c 2 (h + lf 



R 



p(h + iy 



2h(2h+l) (6h(2h + l))* 



R~ 



(17) 



Equating (3) with (17) gives 

2R 2 f" (R) -(h+ l)Rf(R) + (2h + l)f(R) 



+ 



c 2 (h + l) 2 D 2(5{2h +l)(h+l) a « 

R a R 2 



h 



(Qh(2h+1))' 



= 0. 



(18) 
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Solving Eq. (18) yields the PLECHDE /( J R)-gravity as 

c\h+l) 2 R 



.ran = a . /."" + a /."" h , 

A/3(2h + + I)«i2f 



(/i(a - 4) - (a - l)(a - 2))(6/i(2/i + 1))* 
where 

3 + /i ± Vh 2 -10h + l , OAN 
m± = . (20) 

Also A± are the integration constants which can be found from the following boundary 
conditions [9] 

f(R ) = -2R , (21) 
f(Ro) ~ 0, (22) 
where i? ~ (10~ 33 eV) 2 is the current curvature. 

Applying the above boundary conditions to the solution (19) one can obtain A± as 



Ah 



(to± -m T )Ro 



2(3(a - 2m T )(2h + l)(h+ l) a i? f 



(h(a - 4) - (a - l)(a - 2))(6h{2h + 1)) » 



(23) 



Replacing Eq. (19) into (6) and using (8) and (16) one can get the EoS parameter of the 
PLECHDE /(/^-gravity model as 



aX 



2X-1 



(24) 



where 

/?(fe+l)°- 2 , g m+l)- 2 / 1 

6c 2 V 6 /i(2/i+iy " 6c 2 Vi + J ' 1 J 

and z = - — 1 is the redshift. We take ao = 1 for the present value of the scale factor. 

For a < 0, j3 < and — 1 < h, (h ^ 0) one can find that the EoS parameter (24) can 
cross the phantom-divide line. 

For the second class of scale factors, the resulting f(R) is 

f(R) = \+R m+ + \-R m - - 1 ro J + 

AT 



(/i(a - 4) + (a - l)(a - 2))(6/i(2/i - 1)) 

where 



3 - h ± \/h 2 + lO/i + 1 
m± = . (27) 
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3 PLECNA /(fl)-gravity model 

The new agegraphic dark energy (NADE) density is given by [10] 



Pa 



where rj is the conformal time defined as 



V 



3n 2 M 2 



a dt 



(28) 



Jo Ha 2 

The power-law entropy corrected version of the NADE density takes the form [8 

For the first class of scale factors (7) and using Eq. (8), the conformal time rj yields 

rf«df (t s -t) h+1 1 /6/i(2/i + 1)n*J! 



(29) 



r/ 



a (/i + l) ao(/i + l) 



(30) 



(31) 



Note that due to having rj > then — 1 < h. 
Substituting Eq. (31) into (30) one can get 



PA = Ml 



( 3a 2 n 2 (h + 1) 



R 



h+l 



pa 2 (h + l) 



2a 



\(Qh(2h + l)) h+1 
Equating (3) with (32) gives 

2R 2 f"(R) -(h+ l)Rf(R) + (2h + l)f(R) + 



a(fe+l) 

(6h(2h + I))—*" 



Rt 



(32) 



/ aln 2 (h + l) 



Rh+ i_ 



2pal(2h + l)(h + l) 



2a 



a(h + l) 

(6h(2h + 



\h(Qh(2h + l)) h 
Solving Eq. (33) yields the PLECADE /(/^-gravity as 

a 2 n 2 (h + l) 2 R h+1 



o(h+i) 

R 2 



(33) 



f(R) = \ + R r 



\_R m - - 



h 2 (h + 2)(Qh(2h+l))' 



+ 



4f3a 2 (2h + l)(h + l) a R* 



(6h(2h + l))—5—((2 -3a + a 2 ) + 2h{2 -2a + a 2 ) + a(a - l)h 2 ) 



(34) 



where m± are given by (20). Also A± are determined from the boundary conditions (21) 
and (22) as 



1 



(m ± -m T )iC 



2m T R + 



a 2 n 2 {h + l) 2 {h + 1 - m T )R h Q +1 
h 2 {h + 2){6h{2h+l)) h 



2(3a 2 {2h + l)(h+ l) 2a (a(h + 1) - 2m T )R 



a(h + l) 



(2(2/i + 1) - (h + \){h + 3)a + (h + l) 2 a 2 )(Qh(2h + 1))" 
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(35) 



Replacing Eq. (34) into (6) and using (8) and (31) one can get the EoS parameter of the 
PLECNADE /(^-gravity model as 



u - 1 2{h+1] 
oo R — — i — 



aX- 1 
2X-1 



3h 
where 



(36) 



6n 2 ^6/i(2/i+l) 



W2/1+IV 6n 2 h + zl ' 1 ; 



Equation (36) shows that for a < 0, /3 < and —1 < h, (h ^ 0) crossing the phantom- 
divide line can occur. 

For the second class of scale factors we find 



/i 2 (/i-2) 

W/i-iKft-i) 2 ^ 



(6/i(2/i - l))^ i ((2 -3a + a 2 ) - 2h(2 -2a + a 2 ) + a(a - l)h 2 ) 



(38) 



where m± are given by Eq. (27). 



4 /(-R) reconstruction in de Sitter space 

The scale factor in de sitter space is defined as 

a(t) = a e Ht , H = constant, (39) 

which can describe the early-time inflation of the universe [4]. Using Eqs. (5) and (39) 
one can obtain 

Then Eqs. (3) and (4) take the forms 

k 2 p R = -lf(R) + 3H 2 f(R), 
k 2 PR = \f{R)-3H 2 f\R). 

Also the EoS parameter yields u R = ^ = — 1 which behaves like the cosmo logical con- 
stant. 

4.1 PLECHDE model 

For the scale factor (39), using Eq. (40) the future event horizon R h reduces to 



6 



Substituting Eq. (42) into PLECHDE density (14) yields 

R 
12 



Rf'(R) - 2f(R) - c 2 R + 4/3(-|) f = 0. (43) 



This gives 

f(R) = \R*-c*R-^Q f . (44) 

where A is an integration constant. Also from cur = — = — 1 and continuity equation for 
PLECHDE one can get 

Finally one can rewrite (44) as 

r 2 (n - 2) 2 

f(R) = XR 2 - ) > R. (46) 
a[a — 4) 

Note that the term R 2 confirms that the model (46) satisfies the inflation condition [11]. 
4.2 PLECADE model 

For PLECADE, the conformal time r) for the scale factor (39) yields 

f* dt 

Here to obtain rj = rj(R) one cannot replace t by R in (47). Therefore for the scale 
factor (39) one cannot obtain the /(i?)-gravity models corresponding to the PLECNADE 
density (30). To avoid of this problem we set t — > 00 for the upper limit of the integral 
(47). Hence the result yields 

f°° ^ , „ s 1 / alR\~ 112 , t . 



Substituting Eq. (48) into PLECNADE density (30) and using Eq. (41) one can obtain 

a 2 R 
12 



Rf'(R) - 2f(R) - aln 2 R + 4/?(^) f = 0. (49) 



Solving the above differential equation yields 

f{R) = XR 2 -aln 2 R-^ l {^)\ (50) 

where A is an integration constant. Also from ojr = — = — 1 and continuity equation for 
PLECNADE one can get 

a 2 n 2 R ( 12 xf 

fi = -tor($R) ' (51) 

and we finally find 

f { R) = XR'- a ^-f R . (52) 
a(a — 4) 

Here like the model (46) the inflation condition is respected because of the term R 2 . 
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5 Conclusions 



Here, we investigated the power-law entropy corrected versions of the HDE and NADE in 
the framework of f(R) modified gravity. We reconstructed different /(-R)-gravity models 
corresponding to the PLECHDE and PLECNADE models. We obtained the EoS param- 
eters of the PLECH and PLECNA /(^-gravity models. Our calculations show that for 
the selected models, the transition from the quintessence state to the phantom regime 
can occur which is compatible with the recent observations. Furthermore, we studied the 
PLECH and PLECNA /(i?)-gravity models in de Sitter space. We concluded that these 
models can also satisfy the inflation condition. 
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